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I . INTRODUCTION 


This report presents the results obtained for the first 
three tasks under a program whose objective is to analyze 
aspects of the broadband noise generated by rotor or stator 
blades subject to inflow turbulence, in an annular duct. 

The first technical section (Section II) introduces dissi- 
pation, in the form of a generic eddy viscosity, into the 
Green's function relating radiated pressure to blade fluc- 
tuation forces. This adds an element of physical realism 
to the analysis and serves to eliminate the singularities 
or resonance effects which would otherwise occur in the 
output power spectrum. These singularities occur at the 
cutoff frequencies and become dense at higher modes, thereby 
completely dominating the output power spectral behavior. 

The introduction of dissipation eliminates these singular- 
ities and smears out the neighboring peaks so that the 
resulting output power spectrum should be well-behaved. 

The next section (III) then relates the aforementioned 
eddy viscosity to actual dissipation losses experienced by 
waves propagating in a turbulent medium. This is carried 
out in two ways; first by means of a perturbation technique 
of J. Keller for waves propagating in a medium with ran- 
domly variable index of refraction, and second by means of 
scattering losses in a single-scattering treatment of waves 
in a random medium. The second method is carried out by 
integrating over-all scattering angles the single-scattering 
cross-section corresponding first to fluctuations in density 
and sound speed in the background medium and, next, to 
velocity fluctuations in the background medium. These 
scattering cross sections are due to the work of Batchelor' 2 ' 
The results from the two appraoches are shown to agree in 
the scalar-fluctuating case where the background fluctuations 
are the same. They give a wave-number dependence to the eddy 
viscosity, involving the energy spectrum of the background 
fluctuations at all wave numbers less than twice the incident 
wave number. The method of incorporating these results into 
the Green's function for radiated pressure is developed at 
the end of Section III. 

Section IV develops the space-time correlation function 
for the normal component of the turbulent inflow velocity 
fluctuations, assuming an isotropic, frozen turbulent inflow 
velocity correlation tensor. The normal-velocity correlation 

is derived for distinct times and distinct points on 

a single blade or on two different blades. The derivation is 
performed admitting radially varying blade twist angle. 



Then the last two sections V and VI are aimed at develop- 
ing an insight into the qualitative behavior of the blade 
force fluctuation correlations in terms of the normal velocity 
correlation developed in Section IV. Section V derives 
simplified approximate expressions for R^j. in terms of three 
components of the distance function in the case of zero or 
small radial variation of blade twist angle. Then Section VI 
utilizes these results to infer the temporal stationarity 
of the blade-force fluctuations, the homogeneity with respect 
to blade number, and the inhomogeneity with respect to radius. 
These conclusions should then serve as the basis for subse- 
quent developments of output power spectra for radiated 
pressure . 
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II. INTRODUCTION OF DISSIPATION INTO ANALYSIS OF ACOUSTIC 
PRESSURE FIELD 


We consider the linearized (acoustic) compressible flow 
equations for an annular domain in the presence of a uniform 

flow U Q in the axial (x) direction, a force distribution 
f i (per unit volume), and a dissipation in the form of a 
simple eddy viscosity V T . Later the eddy viscosity will 
be considered in detail and made wave-number dependent. We 
retain the simple density-pressure relation 


0(7 - 

despite the presence of dissipation. 


( 1 ) 


The linearized equations are then as follows: 


D; t- (2) 

f, Du> t- >J\ - f, V r AU‘ +• f- ( 3 ) 

where 0 ~ +■ tAo (4) 

bt a* 

Then, eliminating u r and , we find the following 
equation for the pressure perturbation p: 
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We then introduce the cylinder functions 


(6) 

which are linear combinations of and Y ra ( p, m qr) 

such that the normal (radial) derivative vanishes at both 
inner (r^ n ) and outer (r ou ^) radii of the annulus. Further, 
we normalize such that 


•w- 




V" R*ie^, (v) 


(7) 


The R m q then constitute an orthonormal system with 
respect to weighting function r, over the domain (r in , 

r out>* 

Next we operate on eq. (5) with the operator. 
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and define p ( > y c 0) as 
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( 9 ) 


The result of this operation is as follows: 
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Now we recognize that f has only axial and torque 
components f x , f d , respectively, and we perform the inte- 
gration on the right hand side of eq. (10) by parts with 
respect to x Q and 0 0 . The right hand side then becomes 


il 


v„ d~ 


Q fy. ) 


| £ X f-„ v„ 4- ; m Q 




(xc.v'o e 0 fc„ ) (li) 

' / / / 


With the geometry as shown in Fig. 1 and with light loading, 
the resultant force f is closely normal to the blade. Hence 


£.60 

= fcos£fv) (12) 
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where C ( r ) is the local blade twist angle. Then the 
bracketed term in eq. (11) takes the form 


{ \ v; /tv. ) ( l A C £(*v ) 1 M cos £ ( vj) 

(13) 


In the absence of the dissipation term, we can write the 
left-hand side of eq. (10) in the form 


-(i-aC) (14) 

where \, X L are the roots of 



Then p would result from the four-fold inversion of expression 
(11) divided by (14). That is, in the nondissipative case, 
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Now if the integration with respect to A is carried 
out by residues, it is clear that the difference forms 
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will occur in the denominator for each fiu>,m,q. 
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For eacb combination m,q the value Ct) of a) 
which annihilates the radicand in (17) is the well-known 
cut-off frequency which separates propagation from non- 
propagation of the m,q mode. These zeros CO become 
extremely dense as m and q get large. Furthermore, 
when we proceed to form power spectra of pressure p, the 
cross terms (involving are iategrable with 

respect to cj , but the self-terms are not integrable. 

Instead, they lead to infinities or resonances in the output 
or response power spectrum. These resonances become so dense 
at large m and q that they completely dominate any output 
power spectral considerations. 

It is for this reason that the dissipation ( V T ) 
term is introduced. With the dissipation included, we expect 
that these resonances, which occur at the cut-off frequencies 
CJ m q, will be rounded off and, for large m,q smeared to- 
gether, leading to a finite, continuous output power spectrum. 
Moreover, since the origin of the blade fluctuation forces 
is presumed to arise from the turbulent character of the 
inflow to the blade row, therefore it is consistent to assume 
that the dissipation should be of a turbulent character, 
and V-r is to be considered in the sense of an eddy viscosity. 
Subsequently, )) T will be related directly to the dissipa- 
tion experienced by acoustic waves in a turbulent medium and 
will be given wave-number dependence. For the present, 
however, we treat V T as a parameter and return to the full 
expression (eq. 10) for evaluation of the poles in the A 
plane in the presence of nonzero V T 

Now the full equation for the roots of A is as follows: 

(i-m.Ma* +- +- w,<J X 

7 y ciS 

(18) 
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As it stands this is a cubic in X . However, we make use 
of the fact that the dissipation is small and that the 
only places where IV plays any significant role are at 
the cut-off points (or singularities in the nondissipative 
spectrum). Thus we solve eq. (18) for two roots 
as follows: 

ffe treat the V T term as a small correction to the 
quadratic equation for X and lump it in with the constant 
terms, giving A (where ever it appears in this term) 
its nondissipative cut-off value, since these are the only 
points where this term matters. This may be formalized 
somewhat as follows: Consider the generic equation for A 


a,\ u 4 - lo x *-c V T f (\) =o (19) 

where f( A) is any regular nonlinear function of \ . 
Assume is small and expand A in powers of V T . 

(Actually, this should all be done in dimensionless form, 
but the argument would proceed in the same fashion.) 




( 20 ) 


Then the nondissipative roots are 
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And the first order (in Vj- 


) corrections are given by 
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( 23 ) 

We see that the Vf term enters just as it would if it were 
lumped with the c term and given the argument ( 0 ). ‘Now 
we simply give A( 0 ) its values at cut-off inside the function 

f( A (o))> since these are the only points where the Vr term 
enters significantly. 


10 



( 24 ) 



tfhen we apply this recipe to eq. (18) we find 



( ' - Mo 2- ) (25) 


where 
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becomes 


Thus our final expression for the roots X, X, 

' i 



(30) 


Since <\ are always complex, one of them will be in the 
upper half plane and the other will be in the lower half 
plane. In order to avoid confusion we simply call the 
upper root X up and the lower, X lo . 

Note that the form (29) makes it easy to locate X, 
and Ai_ in the upper or lower half planes, depending upon 
the value of U) . In fact, this form, in the limit V r --» o* 
provides the best way (from a physical standpoint) to 
evaluate the pole-locations with respect to the real A axis 
in any residue calculations for the nondissipat ive case. 

Then we may write the expression for the radiated 
pressure field (See eq. 10) in the form 
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Here, the notation 

" oO '' 

\ cl A ^ H (x 'Xu') S fx ) f- H^*x) S f A~A| U )^ 

\ - - to 


is used to signify that all terms involving \ are evaluated 


at 

| A - X ^ \ 

, when 

f X > \ 

respectively . 


' A - Ai„ / 


V X < *o / 



If we concentrate the blade forces onto rotating lift- 
ing lines located in the plane x Q = 0, then 
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Here X U p and A are simply computed from the Aj , 
of eq. (29), depending upon which lies in the (upper, lower) 
half plane of complex X . 
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III. COMPUTATION OF WAVE-NUMBER DEPENDENCE OF DISSIPATION 
OR EDDY VISCOSIT? 

If we consider plane waves propagating at frequency co 
in a quiescent (U Q = 0) free field, in the presence of a 
simple eddy viscosity Vf , we find for an ’'incident'* 
wave j> 0 (referring to eq. (5) with U o =0, f = 0) 


(-k** <i± l t- A k i -to ) ? 0 ~ 
a. 1 J 
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(34) 


or 
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(35) 


= k Lh 

ia 0 


where 


j z < 

k 0 - u > 


a< 


(36) 


Hence the wave amplitude exhibits an exponential attenuation 
as it propagates in the direction of positive x. This is 
of the form, 



Now J. Keller^ considers a wave propagating in a 
random medium whose scalar index of refraction is given by 
the expression 


( 



(38) 


Here fJL is a centered random isotropic process with /A con- 
sidered small relative to unity, and only weakly time -dependent 
relative to the time dependence (d3t) in the "incident signal." 
He finds that the medium can be characterized in terms of 
an effective wave number k which is complex. Moreover, the 
imaginary part of k is shown to be given by 

L - (hIK <m') 1 L'c^ (k ‘) c(k ' 
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where 
and 0 


<M 2 ) 


for 
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is 



is the mean square value or intensity of JL[ 
its power spectrum, and isotropy is assumed 
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e d Sv 
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(40) 


Translating Keller* s k;^ into an equivalent eddy viscosity, 
via eq. (36) leads to 
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Thus we find an equivalent eddy viscosity which de- 
pends upon the power spectrum of the random index of refrac- 
tion at all wave numbers less than twice that of the undis- 
torted incident wave. 

Now G. K. Batchelor^) considers single or Born 
scattering of an incident plane wave in a medium with ran- 
dom fluctuations in density and sound speed. If we temp- 
orarily supress the actual relationship between n 

Q? and d a and consider only the random fluctuation jr“" 
a o 

in sound speed, Batchelor's result can be interpreted as follows 
For a stochastic system satisfying the equation, 
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a centered, random process, an incident 

gives rise to a Born or single scattering 
, in the direction ^ , equal to 
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Here 
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-f is a unit vector, 



is the power spectrum of the fluctuations 
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and 40 is the scattered energy flux per unit volume per 
unit solid angle at direction ^ , normalized by the incident 

energy flux per unit area. 


Thus, since only the if * < term in Keller's index of re 
fraction (38) contributed to k im (39), we should find a cl 
related expression derivable from (£) (43) with Batchelor's 

$ %Cx/e . playing the same role as Keller's $ 


The relation is derived as follows: The integral over 

the unit sphere of clT ) 

v A } 


T = | [ c r(t)it 


(44) 


gives the energy scattered out in all directions, per unit 
volume, per unit incident energy flux. Thus we have an energy 
transport equation for the energy flux S . 
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Therefore we should find 
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where we now assume; as did Keller, that ^ / is isotropic 

and thus depends only on / k Q - k Q ^ f= 2k Q sin ^ . Here 

is the angle between incident and scattering directions, 
k Q and / , 

Integrating (47) with respect to Cp and letting 
2k_sin ® = k* , we find 
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And this agrees exactly with Keller’s result 



playing the role of Keller’s 



(39), 


with 


Now Batchelor ^ shows the relation between fluctuations 
in sound speed and those in density , 


iSii = - It 

Qo A 


(50) 


and gives a more complete form of the stochastic wave equation 

(42) which includes both effects simultaneously. The net 
result, after invoking single scattering, is to modify the 

. -op 

scattering cross section <T (/ ) from the relation given by 

(43) to the following: 
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The angle & is, as before, the angle between k 0 and / . The 
spectrum is now that of the relative density variation. 

Repeating the preceding development for this new scattering 
cross section and again invoking isotropy, we find 
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Hence, finally, the equivalent eddy viscosity for isotropic 
fluctuations in density, with related fluctuations in sound 
speed, is given by 
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( 53 ) 



And again, in this more complete case, the effective eddy 
viscosity depends on the power spectrum of ^f/f for all 
wave numbers less than twice that (k Q ) of the incident wave. 

The relation (53) for a wave-number dependent equivalent 
eddy viscosity holds for fluctuations in the related scalars, 
density and sound speed, in the background turbulent medium. 

A similar, but more complex treatment can be given for the 
single scattering approximation to. fluctuation-velocity- 
induced dissipation. Batchelor^/ also treats this case and 
arrives at a scattering cross section 


<T 



given by (in our notation) 


c r(i ) - z rr <V > k 
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where the notation is as before, with the additional quantities 



mean square fluctuation velocity (in one direc- 
tion) 

fluctuation velocity power spectral tensor 


—9 

We integrate (T over the unit sphere (over X ) as pre- 
viously, to find X for an isotropic velocity fluctuation field 
in the following sequence of steps: 
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Then using 
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r'lq' 



we find, for the equivalent eddy viscosity resulting from 
velocity-induced (single) scattering losses, 



O 


(59) 


Thus, once again, we have an equivalent eddy viscosity 
dependent upon the energy spectrum of the background fluctuations 
(in velocity, for the present case) for all wave numbers up to 
twice the incident wave number. 

The actual eddy viscosity to be used in our expressions 
(e.g. eq. 29) for A should include both contributions (59) 
and (53). If coupling exists between scalar Sf/p^and vector 
Sit fluctuations, then there should also be cross spectral terms 
in Vt . However under isotropic assumptions if we invoke 

an additional assumption of nearly solenoidal velocity fluctua- 
tions, then velocity fluctuations are uncorrelated with scalar 
fluctuations 


<6jSu> - o 


(60) 


We invoke this approximation herein and merely add the contribu- 
tions (53) and(59) to ^ . 

The foregoing has been developed for a quiescent free field. 
To find the value of k Q appropriate for our application (eq. 

29) we proceed as follows: First in fluid-fixed coordinates 
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( S j T') frequency is related to that cD in our 

(casing-fixed) (xyzt) coordinate system by the equations, 

X- ^ fUot 

y ~ i 

f = r 


(61) 


Therefore we have the equivalence, 


£ (Ax-t-o-b) fYx-C +• AU a c -cw>C 

e ~ e * 


(62) 


Hence the relation between and lO is simply, 


and 




(63) 



a u 


(64) 


At cut-off 



/ - M 


L 

c 


( 65 ) 
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and 


i *- M 0 L \ 

/"V 1 

( 66 ) 



This, or rather its absolute value, is the value of k Q to be 
used in applying eqs. (53) and (59) to the evaluation of 

A or > up and h (eq. 29). 

This result can be shown in a different manner. Our wave 
functions are of the form 




satisfy the equation 


~ /^ M <j f ~ 


( 71 ) 


Hence each of the modes, from which our pressure expression 
is synthesized, has a wave number k 0 satisfying 



(72) 


But, from the first of eqs. (15) we have, (in the nondissipative 
case) 


fU, (- X 1 - - M. A +-j£ -Z-M.XoJ 

T a, 1 . 


( 73 ) 


This agrees with our expression (64) for the value k Q to be 
used in evaluating )> T . 
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IV. CORRELATION R xi _ FOR BLADE -NORMAL VELOCITY FLUCTUATIONS 

In forming output power spectra or related quantities 
from the expressions developed in sections II and III, it will 
be necessary to have available the space-time correlation of 
the forces 

on blades K , at radii ri,r 2 and at time instants 

respectively. In order to perform an actual computation 
of such correlations it would be necessary to invert the inte- 
gral equation for three-dimensional, nonsteady compressible 
cascade aerodynamics under sinusoidal gust entry conditions. 

Since this is beyond the present state of the art, we proceed 
as follows: First, in the present section, we compute the 

space-time correlation Rj.i for blade-normal velocity compo- 
nents on one blade or on two distinct blades. Then, since 
the blade force fluctuations are induced by the inflow normal 
velocity fluctuations, we use the results for R u , or simplified 
approximations thereto, to infer certain facts about the ex- 
pected general behavior of the required blade force correlations. 

We refer to the blade-row geometry as illustrated in 
Fig. 1, and develop the correlation R a x in terms of the 
velocity correlation tensor R^j of the inflow. Moreover, we 
assume that the inflow correlation tensor is isotropic in fluid- 
fixed coordinates and frozen in casing -fixed and blade-fixed 
coordinates. This means that, in fluid-fixed coordinates 

(M-I/T) , R i;) is given by <3 > 


R; j = +• 3'fJfci 

r 

with f - 

and f'~ f ~ ( t V< 1 ‘ 1) ) 



R i - is considered normalized, so that f and g are dimensionless, 
with unit values at 0. 


f (o) - ° ) " 1 


( 75 ) 


Now 

x - t; x U«r 

y 

*- - i 

ir ~ L (76) 

and thus for the V th blade lifting line located in the 
plane x = 0, 

V~ a ° b 

x = t 

( 77 ) 

Therefore, for two blades ^ and two different times 
ti,t 2 and radii ri,r 2 , we find the following: 



- U-tj 


( 78 ) 
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where 



Also 

f ~ | ?| - UM/ I" V *, 1 h - U, VLCoS A 


where 


( 79 ) 


( 80 ) 


A-ALfr a + *TT ^ (81) 

N 

From Fig. 1 the normal fluctuation velocity component u^ 
can be expressed as follows: 

U L ~U A 5i>i^ ■*- \Aj_ Co5 £ CoS (-Q- f +- Itr l//rt) 

- (-air +- ud//n) 

_ U, S W\Z +- Co5e[ UiCos(^ll'+ WVM ) 

- UiSl \i\(-Slir f- Zf TV/ii ) ] 

( 82 ) 
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where 


e - z(v) = 


£ -Lev) = blade local twist angle as 

defined by reference to Fig. 1. 

Since we are going to form a second order correlation and since 
( u x>Uy»u z ) have zero mean values, therefore, it is unnecessary 
to correct for mean velocities -*1 t and U Q . 

Now we form R^ x from eq. (82) evaluated at (r-j^t^) on 
blade and at on blade V 2 . 

< V ,l > R xL * < V' z ? | |?„ £, Sm i- case, C oS € 2 ( c, c l 

f Sw £, C oS£i Ci' 


H Sh CoS £, ^ c, - |?^ 5 , J | 


(83) 

Here the are given by (74) with f and f as given by (78) 
and (80), respectively. < is the (unidirectional) mean 
square velocity fluctuation. 
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Thus 


R„ ~ ( ^ ^(f) 

s x 

R 33 ~ ( )- £ j(f)) l V.s, ~V\ S if h ^f) 

y L 


R«*. = = ( (-ff) '^f) )(~^t,i)(v,C l -\r tr C L J 

r 

^•S = ^3< ~ ( ^f)~ ^(f) ) (~^k<)(V, S,-V\Sz.) 

- r ~ 

f x 

(84) 



( 85 ) 


where f and A are given in eqs. (80) and (81). For zero 
or small twist variation 


£( *-) - £, = = t- 


( 86 ) 


and 


K 


i-1 


( .. i 

^ j £ t- Co 


sA a oS Z t ) 


(f ~ \ ) f S*.^£ ^ j-cos ? £ V\ 


— Sl^ £ CoS £ 


t*l (V, Hr L ) j 


(87) 
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( 88 ) 


we may rewrite eq. (80) for y 2 in the form 

u/t/ f z\r,v-^(i~c*s A) 

-U C L ^ t~ +- 4- (v M z - )ji *\ A 

«" *■ (89) 

If we assume that f(f> ) g( f ) decay sufficiently fast with 

P that we may neglect terms in p 2 of higher than second order 
in difference quantities ( )d, then we may approximate j 1 *" by 
the expression, 

O 1 - ^ f i- f 



Iso the constant- £.form (87) for R iJL 
by 


can be approximated 


R aj . ^ +■ ^ sWe Ujfy 1 +W<e 

- 2_ $ £_coS £ U^jr^YVn") 


m£ 4 - ^ c 05 £ 

$ 


2 . 




y?) J- y^)-'s(f) (-LL^S'^f i-cose r ^j 

" f‘ 


Z. 


(91) 





as follows: 



(v, " r i) L ~ (radial) 

^J. — (‘ Mo'htS fv* £ 4. t c & S'£. ) 

u 

•=" (*“ Wol^'SK c f r^dCcSe 

/ 

5 (, ~ ( ~ . c<?5£ ~ t^S-1^4 $<>1 £ j 

- ( “ C 6 $ £ 


(blade-normal) 

(92) 

(cbordwise) 
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Then 2 , as given by approximation (92) is equal to 


(93) 

and 

R lL ~ Ckffl f- (-f-ff) -<$(?)) Sfc* 

J f Z 

- Qi Wlti 1 *- O *■ f(?) jkL 

^ f * 

(94) 

For solenoidal velocity fluctuations, ^ 2 ) 

i 

% 

Hence, in this case 

R^s Wf) 

= f 1 f 1 

( 96 ) 


[ ’ft) f f-fjJ 

05) 

+• 1 

*- f/) 

■% O ' 
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For example, if f (j> ) is approximately expressible by 


fr f ) - e f/A 

(97) 


and if the condition (95) is met, 
Then 



f..W ) 

2j>A 


(98) 


where A is the integral scale for f. The correlation 
6L leads to a k -2 type power law in the power spectrum 

E(k) and hence is not a bad approximation in the universal 
equilibrium range where k“ 5//3 is the correct dependence. The 
simple exponential fails at y = 0 where it has a discontinuous 
slope. 
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VI. INFERRED BEHAVIOR OF BLADE FORCE CORRELATIONS 

In accordance with the plan outlined at the start of 
Section IV, we use the simplified forms of Rxx developed in 
Section V to infer something about the general behavior of the 
blade-force correlation in the absence of a full inversion of 
the integral equation for three-dimensional sinusoidal gust 
entry for a cascade in subsonic compressible flow. 

We define the normalized correlation Y' as follows 

'i' 

(99) 

The forms given in eq. (94), together with the definitions 
(92) of ?i_ and , are the most suggestive for purposes 

of inferring the general behavior of Y' . First Y' should 
obviously depend on r^ and r 2 through both r m and r^, since 
and P„ exhibit this behavior. Second, Y' should depend on 
tj and t 2 only through t d , since t m enters nowhere in f 
or Ru and since the boundary value problem for the nonsteady 
blade pressure distribution is invariant under translation in 
time. Thus Y' is stationary (independent of t m ) . Finally, Y' 
should depend on blade numbers V t and y x only through the 
difference ]/^ , since only (and not ) enters into 

R i (through A ) and since the boundary value problem for 
pressure is also invariant under translation in 9 . Hence 

we expect 

= ^ , a, (- ) ) 

( 100 ) 

and A , in turn, depends on ^ and . Thus Y' should be 
stationary in time and homogeneous in blade number, but not in 
radius.- 
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In order to invoke the best available estimate of 
a wave-number dependent transfer function between upwash 
and blade lift per unit span, we revert temporarily to 
a two-dimensional cascade geometry and utilize the single- 
wing transfer function developed by Filotas^^ for airfoil 
response in a gusty atmosphere, together with the inter- 
blade cascade effects developed by Lane Friedman' 5 ' for 
subsonic compressible oscillatory aerodynamics. The 
Filotas analysis treats three dimensional sinusoidal gust 
entry in terms of the effective inflow wave number ampli- 
tude and direction. The resulting transfer function 
relates lift per unit span to upwash as follows: 

f # b. 0U(o,o,O,f)T 


Lift 

unit span 


Here 'V(o,o/ 0 ,t) is the inf low upwash at midchord at time t. 
Then T is given as 





TTpfu-^CoSP) 

1 V- ITT V, b,(l+*£Co5(i) 


(i v- rrkk ()*- sm‘ l p ttU, b*cos ($)j 


( 101 ) 


Here k Q is the magnitude of the blade-plane component of 
inflow turbulent (Fourier component) wave -number vector 
k Q and p is the angle between k Q and the blade mid-chord 
line. If we express the turbulent inflow velocity 
(V(X) in the usual Fourier-Stielt jes form 
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(102) 


u*Yx) = 




then for each Fourier component , dzW , k Q is given by 

V 0 s l< - H ( hJ * U ) 

( 103 ) 

. — • 

where N is the unit normal to the blade plane. Then the 
angle is the angle between k 0 and the blade mid-chord 
line. 

The interblade dependence of the transfer function 
may then be estimated by use of the oscillatory aerodynamic 
influence coefficients developed by Lane and Friedman^®) 
using an interblade phase lag angle (r which is expressible 
as 

-A 

0~- tc 'l 


( 104 ) 

Here J? is the vector from the mid-chord point in the n th 
blade to the mid-chord point on the (n+l) th blade, in a 
plane normal to the blade mid-chord lines. 

In this way the best available estimate of wave- 
number dependence and blade-to-blade influence can be 
incorporated into the blade force correlations (or the power 
spectrum thereof) for use in subsequent radiated-f ield 
power spectral estimates. 
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